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Abstract. We give for a plasma with a history of several 
expansion and contraction phases an analytical model of 
the evolution of a contained relativistic electron popula- 
tion under synchrotron, inverse Compton and adiabatic 
energy losses or gains. This is applied to different sce- 
narios for evolution of radio plasma inside the cocoons of 
radio galaxies, after the activity of the central engine has 
ceased. It is demonstrated that fossil radio plasma with 
an age of even up to 2 Gyr can be revived by compression 
in a shock wave of large-scale structure formation, caused 
during the merging events of galaxy clusters, or by the 
accretion onto galaxy clusters. We argue, that this is a 
highly plausible explanation for the observed cluster ra- 
dio relics, which are the regions of diffuse radio emission 
found in clusters of galaxies, without any likely parent ra- 
dio galaxy seen nearby. An implication of this model is 
the existence of a population of diffuse, ultra-steep spec- 
trum, very low frequency radio sources located inside and 
possibly outside of clusters of galaxies, tracing the revival 
of aged fossil radio plasma by the shock waves associated 
with large-scale structure formation. 



Key words: Acceleration of particles - Shock waves - 
Galaxies: intergalactic medium - Galaxies: active - Galax- 
ies: clusters : general - Radio continuum: general 

1. Introduction 

The radio cocoons of radio galaxies become rapidly un- 
detectable after the central engine of the active galactic 
nucleus ceases to inject fresh radio plasma, or even earlier 



(Kaiser et ah, 2000, and references therein). Radio plasma 
consists of relativistic electrons and magnetic fields, which 
are responsible for the synchrotron radio emission, and 
possibly relativistic protons and/or a non-relativistic ther- 
mal gas component. 

Although undetectable, aged radio plasma should still 
be an important component of the inter-galactic medium 
(IGM). We investigate the possibility of reviving patches 
of such fossil radio plasma (also called radio ghosts, Enfilin 
(|l999a|) ) by compression in a shock wave produced in the 



IGM by the flows of cosmological large-scale structure for- 
mation. Such shock waves re-energize the electron popu- 
lation in the fossil radio cocoons, which can lead to ob- 
servable synchrotron emission. This is probably the mech- 
anism responsible for the, so called, cluster radio relics, 
which are patches of diffuse, sometimes polarized radio 
emission typically found at peripheral locations in clusters 
of galaxies. These cluster radio relics can not be simply 
relic radio galaxies, as their name suggests. The spectral 
ages of the electron population are usually too short to 
admit even the nearest galaxy to have been the parent ra- 
dio galaxy, which has moved to its present location with 
a velocity typical for cluster galaxies. Therefore, a recent 
enhancement of the nonthermal radio output of the cluster 
relic sources is mandatory. 

For reviews on diffuse cluster radio emission (cluster 
radio relic and halos, which arc believed to be distinct 



phenomena), s ee Ja ffc ( 1992 ), Fcrctti & Giova nnini (1996 ), 
Enfilin ct al. (|l998|), Fe retti (|l99S| ), Enfilin fll999b|) , and 
Giovannini et al. ( 1999 ). 

The connection between the presence of a shock wave 
and the appearance of the cluster radio relic phenom- 
ena was assumed to be due to Fermi-I shock acceleration 
of electrons, in Enfilin et al. (1998) and Roettiger et al. 
( |1999| ). While this process might work within the normal 



IGM, several arguments favor the possibility that indeed 
old fossil radio plasma is revived in the case of cluster relic 
sources: 

— Cluster radio relics are extremely rare, whereas shock 
waves should be very common within clusters of galax- 
ies. The dual requirement of a shock wave and fossil 
radio plasma, for producing a cluster radio relic, would 
be an attractive explanation for the rareness of the 
relics. 

— Fossil radio plasma with existing relativistic electron 
population and fairly strong magnetic field appears to 
have ideal properties to be brightened up during the 
shock's passage. 

— The cluster radio relic 1253+375 near the Coma cluster 
of galaxies appears to be fed with radio plasma by the 
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nearby galaxy NGC 4789 (see Fig [l] and scenario C in 
Sec. |). 

But, if indeed the fossil radio plasma and not the nor- 
mal IGM were to become radio luminous at a shock wave, 
the expected very high sound velocity of that relativistic 
plasma should forbid the shock in the ambient medium to 
penetrate into the radio plasma. Thus, shock acceleration 
is not expected to occur there. Instead, the fossil radio 
plasma would get adiabatically compressed, and the en- 
ergy gain of the electrons is expected to be mainly due 
to adiabatic heating. It is the purpose of this work to 
demonstrate that this process is sufficient to account for 
the cluster radio relics. 

Sec. ^ and || discuss the expected properties of the 
components of this model: shock waves and fossil radio 
plasma. In Sec. ^ we give a mathematical formalism to 
describe the evolution of the radiation spectrum during 
the history of a radio cocoon. The different phases of this 
history are described in Sec. ^ and applied to three differ- 
ent scenarios in Sec. ^. The discussion of the merits and 
potential deficiencies of this model can be found in Sec. [7J. 



NGC 4789 and Coma's cluster radio relic 1253+275 at 608.5 MHz 
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2. Structure Formation Shock Waves 

The flows of the cosmological large-scale structure forma- 
tion are predicted to produce frequently shock waves at 
the boundaries of clusters and filaments of galaxies and 



during cluster merger (Quilis et al., 1998; Miniati et al 



2000b). But only in a few cases cluster merger shock waves 
could be detected directly as temperature structures in the 



Fig. 1. Radio map of the cluster radio relic 1253+275 in 



the Coma cluster from Giovannini et al. (1991) and its 
probable source of radio plasma: the narrow angle tailed 
radio galaxy NGC 4789. The radio tails of NGC 4789 indi- 
cate the direction from which this galaxy came: the cluster 
center which is located outside this map. A radio polariza- 



tion of 27 % was detected from the radio relic (Giovannini 



X-ray Emitting cluster gas. The best examples are in Abell et al., 199l[ ), but not from the tails of NGC 4789 feeding 

the relic with radio plasma. This is a clear signature of the 



2256 flBrjej and Henry, 1994| ) and in Abell 1367 ( ponnelly 
et al., 1998| ). The reason for the small X-ray detectabil- 
ity of shock waves is the necessary simultaneous spatial 
and spectral resolution of the instrument, combined with a 
sufficiently large collecting area to yield adequate photon 
statistics from the often peripheral shock waves in clus- 
ters. For instance one of the best studied cluster, the Coma 
cluster, is proposed to contain a shock wave at the location 
of t he diff use radio source 1253+275 (see Fig. [|; En filin et 



alignment of magnetic fields enhanced in the shock com- 
pression. The polarization vectors are consistent with a 
shock wave oriented para llel to the main axis of the radio 
relic (jEnfilin et al., 199^ ). 



[gicmicniec and Ostrowski, 200C ; Miniati et al., 2000a 



et al. 



gamma ray background production (Loeb and Waxmanj 
al. pp^ihb work). The available X-ray map gjj-j^ and in bi the ent ic mst of thc lG ^\ 



hIiowh an extension of the X-ray emitting gay 

in the direction of this radio source, but the sensitivity is 
too low to reveal any structure at its location. The feasibil- 
ity of shock wave detections with new instruments as the 
XMM experiment is estimated in Tozzi et al. (2000). It is 
interesting to note that at the locations of the two above 
mentioned X-ray detections of shock waves, cluster radio 
relics are observed (Enfilin et al., 1998). This indicates 



( fTozzi et al., 200C| ). 



that these sources are indeed tracers of shock waves. The 
importance of a detailed understanding of these tracers of 
shock waves is underscored by the shock wave's possible 



|et al., 1997 




Ryu et al., 1998), cosmic ray acceleration 


(Norman et al., 1995; 


Kang et al., 1996; 


Kang et al., 1997; 



3. Fossil Radio Plasma 

In this work it is assumed that the radio plasma in the 
lobes of radio galaxies remains mostly intact after re- 
lease, forming the proposed radio fossils or ghosts. This 
is specified and discussed in the following. Important for 
the presented model is that the magnetic fields and ultra- 
relativistic electron population continue to remain within 
the lobe for up to 2 Gyr. Transport or decay processes of 
magnetic fields, and relativistic particles populations need 
to be slower than this time scale. 
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Clearly, the light radio plasma is subject to buoyancy- 
induced motions in the cluster gravitational potential. The 
action of buoyancy is e.g. able to explain the morphol- 



ogy of the eastern radio lobe of M87 ( Churazov et al 
2000|)~|But since this does not separate the electrons and 
fields, the fossils would not get disintegrated. Neverthe- 
less, relative motion at a plasma-field interface (e.g. the 
surface of a buoyant fossil) can produce enhanced mixing 



with the radio plasma, since both phases are still sepa- 
rated by magnetic fields. Further, there is observational 
evidence that the density of entrained gas cannot be high 
in radio lobes and in cluster radio relics. Many lobes and 
relics exhibit radio polarization, which has sometimes val- 
ues close to the upper limit allowed by synchrotron theory 



(Lerche and Parker, 1967) due to finite gyro-radius effects 



of the ions in the embedding gas. This should produce a 
mixing layer of the order of several gyro-radii of the ions 



(Harris et al., 1993; Saripalli et al., 1996, e.g.). No intrin- 
sic Faraday rotation, or intrinsic Faraday depolarization 
is detected in these systems. In the cases of detected de- 
polarization it seems to be consisten with beam depolar- 



izatio n, or the Laing-Garrington effect ( Garrington et al 



(10 ..10 kpc). But this would soften the magnetic 1988 ), which is external to the source. This argues against 



boundary and therefore inhibit the mixing (Lerche and 



Parker , 1967 ). The entrained gas in that boundary layer 



might be able to propagate into the fossil following the 
field lines, but the large sizes of the fossils (10...10 3 kpc) 
would forbid their disintegration due to the gas leakage 
within the considered time scale. The relativistic parti- 
cles contained in the old radio plasma have much larger 
gyro-radii (lO^kpc (E/GeV) (B/fiG)^ 1 ) than the exter- 
nal gas particles. Scattering at magnetic irregularities al- 
lows the gyro-center of the particle motion to be displaced, 
leading to cross field diffusion required for particle escape. 
But even if the combined effect of cross field diffusion and 
rapid parallel diffusion along disordered field lines is con- 
sidered (anomalous diffusion) , only very special turbulent 
conditions allow an efficient escape of relativistic parti- 



cles from radio plasma (Enfilin, 1999a). Strong turbulence 
itself is probably more destructive by shredding the mag- 
netized plasma into smaller pieces. But it has been argued 
flEnfilin, 199"9^ ; |Enfilin and Kaiser, 2000| ) that this should 
stop on scales where the energy density of the turbulent 
energy cascade is insufficient to overcome magnetic forces. 

From purely an observational viewpoint there already 
exists good evidence that clouds of synchrotron plasma 
can survive and retain their identity within the ICM for at 
least up to ~ 10 8 yr, in the form of relic radio sources show- 
ing a relaxed structure but well defined boundaries. These 
fairly robust age estimates are based on interpreting the 
observed radio spectral steepening in terms of the inverse 
Compton losses of the relativisitic electrons against the 
ubiquitous cosmic micro wave background ph otons (e.g., 
Komiss arov fc Gubanov ( |1994| ); Venturi et al. fll998j ); Slee 
& Roy ( |1998D ). 

Finally, we note that other energy loss mechanisms 
of the electrons, as Coulomb and Bremsstrahlung losses, 
are inefficient at the considered energies (1-10 GeV), even 
if a dense gas component were to exist inside the radio 
plasma ( |Enfilin and Kaiser, 2000 ). If indeed such a cold 
gas component were present inside the radio lobe, Parker- 



type instabilities (Parker, 1966) might occur, helping to 
shred the lobe into smaller pieces. Also Rayleigh- Taylor- 
type instabilities of the surface of the light radio plasma 
might occur, if not supressed by the tension of the mag- 
netic field lines. But it is questionable if these instabili- 
ties allow a microscopic mixing of the environmental gas 



a dense thermal plasma component. Nevertheless, we dis- 
cuss briefly in Sect. some possible implications of strong 
microscopic mixing for the proposed model. 

4. The Formalism 

The dimensionless momentum of an ultra-relativistic elec- 
tron p = P c /(m e c) within the radio plasma changes due 
to synchrotron losses proportional to the magnetic energy 
density ub, inverse Compton (IC) losses proportional to 
the cosmic microwave background radiation field energy 
density uc, and adiabatic losses or gains connected to the 
change of the volume V of the radio plasma: 



dp , , 2 1 1 dV 

-=a (u B +u c )p +3^P, 



(1) 



where ao — |ctt/(tooc). We do not consider 
bremsstrahlung and Coulomb losses in view of the very 
low particle density within the radio plasma. We further 
assume sufficient pitch angle scattering to keep the elec- 
tron pitch angle distribution isotropic ( Jaffc and Pcrola 
1973). After introduction of the compression ratio 



Cit) = V /V(t) (2) 

and a temporary change to the variable p(t) = 
C(t) -1 / 3 p(t) Eq. [I] is easily integrable, yielding finally 

Po 



C(t)-i +p /p*{t) 



(3) 



Here, we defined the characteristic momentum p*(t), 
which is given by 



1 



(J) 



If the change in volume can be approximated to be a 
power-law in time, 



vn) = u [£) ,orC(f) (£ 



(5) 



which we will assume in the following, an analytic solu- 
tion to Eq. [I] was already given in Kaiser et al. (1997). 
We assume further the photon energy density to be con- 
stant, which is reasonable in our application since it is 
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dominated by the cosmic microwave background, which 
does not change significantly on the time-scales considered 
here. Further, we assume that the magnetic field energy 
density scales as 



U B (t)=U B ,0 (V/Vq) 



-4/3 _ 



= u B .o 



(*Ao) 



-46/3 



(6) 



as one would expects for an isotropic expansion of the 
magnetized plasma. These assumptions allow one to eval- 
uate the integral in Eq. |]. We exclude the cases 6 = 3 and 
6 = 3/5 for simplicity only (they produce logarithms in- 
stead of power-law relations), and write the characteristic 
electron momentum in a compact form: 



P*(t) 



dot 



1-56/3 



gl/3_Ql/6 

1-6/3 



(7) 



(C = C{t) for brevity). It is obvious that the synchrotron- 
and IC- cooling produces a sharp upper cutoff in the elec- 
tron distribution f(p,t) at p*(t) even if the original dis- 
tribution /o(po) at time to extended to infinity. The elec- 
tron density per volume and momentum f(j>,t) dpdV for 
p < p* (t) is given by 



dp 



where 



Po(p,t) 



P c(ty 



(8) 



(9) 



(10) 



1 -P/P*(t) ' 

If the original distribution function was a power-law 

/o(po) = foP a ° 
for Pmino < Po < Pmaxo the resulting spectrum is 

/( Pl t) = / c(t)^p- ac (i-p/p*wr~ 2 (ii) 

for p min (t) = p(p min0 ,t) < p < p max (t) = p(p ma xO, *)• 

We are interested in the situation where several phases 
of cooling characterized by different expansion or con- 
traction rates and durations shaped the electron distri- 
bution. We write p\ — p(t\) for the momentum of an elec- 
tron originally at p after phase 1, which is characterized 
by the compression during this phase Cqi = C[t\) and 
P*o\ = P*(t\), P2 for the momentum of the same electron 
after phase 2, characterized by C\i and p*i2, and so on. 
It is straightforward to show that the final electron mo- 
mentum p n after n such phases can still be written in the 
form 



Pn{Po) 



where 



(C „) 3 +Po/p*On 



Co n — Y[ Ci-n — Y\ 



Vi-l 

Vi 



(12) 



(13) 



is the compression ratio between original and final config- 
uration and maximal final momentum is given by 



1 



P*0n 



E 



(Coi-i)J 

P*i — 1 i 



(14) 



The effects of the individual cooling phases sum up 
weighted by (Coi) 3 . Thus, whenever the radio plasma is 
most extended, cooling is inefficient. 

It remains to provide the parameters describing the 
different phases. Suppose we want to describe a phase i 
where the expansion or compression is described by &j, 
and two of the following three quantities are given: n, the 
time scale of expansion, Ci—n the compression ratio dur- 
ing the phase, and the duration of the phase. Theses 
quantities are related via 



Ci. 



We get 



(1 + Au/n) 



-bi 



/3_ C l/b 

-5 6/3 



U Bi 



C l/3_ C l/b 

1-6/3 



u c 



(15) 



(16) 



with C = C'i-n and t{ — Ti+Ati for brevity. The magnetic 
energy density at the beginning of phase i is that of the 
end of phase phase i — 1: 



UBi-l 



«B0(C 0i -l) 4/3 . 



(17) 



The resulting electron spectrum from an initial power-law 
distribution is 



fi(p) = f c 0i 3 P - a < (i- P /p* 0i y 



(18) 



for Pmini = Pi(p m in0) < P < Pmaxi = Pi(Pmaxo) and 

fiip) — otherwise. The synchrotron emission at a given 
frequency v is 



L vi = c 3 B i V i I dpfi(p)F(v/vi(p)), 



(19) 



where C3 = v / 3e 3 /(47rm e c 2 ) and the characteristic fre- 
quency is Vi{p) — 3eBip 2 /(4:iTm c c). The dimensionless 
spectral emissivity of a mono-energetic isotropic electron 
distribution in isotropically oriented magnetic fields F {x) 
can be approximated (Enfilin ct al., 1999): 



^) ^/3 (V^ 

{) r(ii/6) p 



, U x 7/8 



(20) 



In reality, after shock passage an originally isotropic en- 
semble of field lines gets partially aligned with the shock 
plane. This is also true for the unshocked radio plasma, 
since its morphology gets significantly flattened during 
compression (see phase 3 in Sec. [5]). This would produce a 
radio polarizat ion and a luminosit y which depends on the 
viewing angle ( Enfilin et al., 1998|) . As long as we are only 
calculating the total luminosity of the radio cocoon/relic, 
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this can be ignored. But, in case one wants to know the ex- 
pected flux, one has to correct for the anisotropic emission 
pattern of the cluster radio relics. Fortunately, the degree 
of radio polarization can be used to determine the viewing 



angle with respect to the shock plane (Laing, 1980; Enfilin 
et al., jggf ). 

The upper cutoff in the electron distribution at p*o n 
produces a cutoff in the synchrotron spectrum near v m = 
v c {p*o n )■ But since F(x) has a broad maximum even a 
sharp cutoff in the electron spectrum gives a soft cutoff in 
the radio, with significant flux above v m . 



5. The Model 

Between the release from a radio galaxy and the reappear- 
ance as a cluster radio relic the radio plasmon undergoes 
several different phases of expansion and contraction: 

Phase 0: Injection. The radio galaxy is active and a 
large expanding volume is being filled with radio plasma. 
The expansion of this cocoon is likely to be supersonic 
(with respect to the outer medium) and therefore bo = 
9/5, if we assume that there is no gas density gradient 



in the vicinity of the radio galaxy (Kaiser and Alexander 



1997 ). We assume that the injection occurred with a time 
constant of r = 0.15 Gyr, which is the typical age of a 



radio source at the end of nuclear activity (Alexander and 
Leahy, 1987, for typical ages). The particle population is 



kept close to a power-law distribution by injection of fresh 
electrons. We assume a spectral index of a c — 2.5 and an 
upper cutoff at p max o = 10 5 . The momentum cut-off pro- 
duces a radio cut-off above 42 GHz (B/fj,G), but we note 
that our results are not very sensitive to the choice of this 
parameter. More realistic electron spectra at the end of 
phase could be constructed by superimposing the spec- 
tra of the electron populations of different ages, as in the 
models of Kaiser et al. (1997). But, for the present pur- 



pose of demonstrating that the fossil radio plasma can be 
revived by compression, our simplified treatment should 
be sufficiently illustrative. 

Phase 1: Expansion. After the central engine of the 
radio galaxy became inactive, the radio cocoon might still 
be strongly over-pressured compared to its gaseous envi- 
ronment (Begelman and Cioffi, 1989). If this is the case 
a Sedov-like expansion phase exists with b\ = 6/5. Mo- 
mentum conservation of the expanding shell around the 
cocoon forces the expansion rates at the end of phase 
and beginning of phase 1 to be the same, leading to 
n = 61/60 To = 2/3 tq. The expansion will significantly 
deviate from 61 =6/5 at the moment when the internal 
pressure drops to a value comparable to the environmen- 
tal pressure. We simplify this behavior by assuming that 
the expansion is Sedov-like until pressure equilibrium is 
reached. The radio cocoon probably becomes undetectable 
during this phase, and therefore becomes a fossil radio co- 
coon or a so called radio ghost. 



Phase 2: Lurking. Pressure equilibrium with the en- 
vironment is reached, and the volume of the radio plasma 
remains more or less constant. Thus C12 ~ 1 and 6 » 0. 
Taking the appropriate limit r 2 — > 00 of Eq. [l^ gives the 
well known result 



P*12 



[a (u B2 + u c ) At 2 ] 



(21) 



Due to the previous adiabatic energy losses of the elec- 
trons, they reside at low energies during phase 2. Their 
radiation losses, which strongly depend on the particle en- 
ergies, are therefore strongly diminished. Additionally, the 
synchrotron losses are further reduced due to the weaker 
magnetic field during the expanded state of the radio co- 
coon. The adiabatic losses are reversible, and will be re- 
versed during the subsequent compression phase, whereas 
the radiative losses are irreversible. Since the latter are 
suppressed during this phase, the radio ghost state can be 
called the energy saving mode of a radio cocoon. 

Phase 3: Flashing. The fossil radio plasma gets 
dragged into a shock wave of cosmic large-scale structure 
formation, e.g. at the boundary of a cluster of galaxies or 
in a galaxy cluster merger event. While its thermal en- 
vironment gets shocked, the radio plasma itself is only 
compressed adiabatically, due to the much higher internal 
sound speed. The electron population and the magnetic 
field gain energy adiabatically, leading to a steep enhance- 
ment of the synchrotron emissivity. 

This phase has a duration of the order of 
At 3 = y 2 1/3 /w s hock2 ~ 0.1.. .lMpc/(300... 3000 km/s) w 
30 Myr...3 Gyr. The pre-shock flow velocity in the shock 
frame, i> s hock2, is related to the pre-shock sound speed 
c S 2 = \Jl Pij{n c (rap + m c )/2) via 



3»«ka = f^ (t-1 + (7 + 1)^ 



(22) 



( Landau and Lifschitz, 1966 ), where 7 = 5/3 is the adi- 
abatic index of the thermal gas. The compression factor 
of the relativistic plasma is high, and can be calculated 
from the assumed pressure jump P3/P2 of the surround- 
ing thermal gas, assuming pressure equilibrium before and 
after the shock passage: 



C23 — (P3/ P2) 



(23) 



In order to guess the parameters 63 and T3 a rough 
picture of the compression process is required. The sound 
speed within the radio plasma should be much higher than 
even in the post-shock thermal environment, (it could be 
up to c/y/3, if the plasma is fully relativistic), so that an 
instantaneous response to environmental changes can be 
assumed. During the shock passage, the cocoon is exposed 
to the high thermal pressure of the post-shock gas on its 
down- and to the ram-pressure of the pre-shock gas on 
its up-stream side. But, on the remaining surface the co- 
coon is only subject to the (much lower) up-stream gas 
pressure. The relativistic plasma will therefore start to 
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1. 





0.595 
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8.5 
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2. 
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22. 
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0.19 
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0.033 
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Phase 2 
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0.0036 
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0.13 


-0.17 


15.9 


2. 


2.83 


0.0021 


2.4 


0.14 


16500. 


2.7 


1.4 


0.12 


Phase 4 


0.1 


CO 


1. 





2.83 


0.0021 


2.4 


0.14 


7860. 


0.62 


1. 


0.02 



Table 1. Parameter of the three szenarios. The rows labled with 'Phase 0' give the initial parameters of the radio 
plasma at the end of the activity of the radio galaxy. The rows below ('Phase 1-4') give the subsequent evolution, due 
to aging and expansion or compression described by the left block of parameters. Details can be found in the text. The 
luminosity units allow an easy conversion to the observable flux: for a source in 100 Mpc distance the flux of 1 Jansky 
corresponds to a luminosity of 1.2 • 10 31 ergs -1 Hz -1 . 



expand orthogonal to the gas flow, producing a flattened 
pancake- like morphology. In order be able to expel the am- 
bient gas sideways additional internal pressure compara- 
ble to the ram pressure of the expelled material is needed. 
This pressure is produced by compression. The process of 
flattening stops when the ram pressure of the swept-up 
material at the expanding edges of the pancake is of the 
order of the ram pressure of the incoming flow. This im- 
plies that the ratio of the diameter to the thickness of the 
expanded cocoon is roughly 4 for a strong shock. 

The compression is slow in the beginning, and rapid to- 
wards the end when the cocoon is significantly flattened. 
We mimic this by setting 63 = 2 and a negative T3, accord- 
ing to Eq. We favor this over the more intuitive choice 
of a positive T3 and large negative 63, since it describes 
the process of compression more realistically. 

Phase 4: Fading. The radio plasma is in pressure 
equilibrium with the post-shock medium, which should 
provide roughly a uniform environment: 64 = 0. The radio 
emission of the relic now fades away due to the heavy 
radiation losses. 

6. The Scenarios 

We follow the evolution of a cocoon of radio plasma for 
three different scenarios (A, B, and C) and calculate the 
resulting radio emission. These scenarios are envisioned 
to illustrate plausible situations, rather than try to re- 
produce precise parameters for some known cluster radio 
relics. Scenario A and B are chosen to represent extreme 



cases: In A, the relic is located at the center of a galaxy 
cluster, while in B, its location is near the cluster bound- 
ary, i.e., in the proximity of the accretion shock wave. In 
both scenarios the duration of phase 2 was chosen to be 
so long that the shocked radio plasma could be barely ob- 
served as a weak ultra-steep spectrum source. Since these 
scenarios are rather extreme, we demonstrate with sce- 
nario C that a shorter phase 2 can result in a moderately 
steepened spectrum of the cluster radio relic. 

For the cases A and B, we assume the initial cocoon (at 
the end of phase 0) to contain magnetic fields, relativistic 
electrons and protons with energies of E B / c / p = 10 60 erg 
each, and i?B/e/p = 10 58 erg in case C. The three compo- 
nents produce a relativistic isotropic pressure of 

_ E c + E p + E B , , 

r cocoon0 — o ' \ ) 

6 Vo 

We assume a spectral index of a c — 2.5 and a rather 
high cutoff in the electron spectrum at p max o = 10 5 , which 
has no significant influence on the conclusions. The lower 
cutoff is p m in o = 10 in scenarios A and B, and p mm o = 100 
in C. The lower cutoffs only affect the normalizations of 
the radio fluxes, not the spectral shapes. 

The different scenarios are described below in detail. 
The parameters for the different scenarios are listed in 
Tab. |] and the resulting radio spectra at the end of the 
different phases are shown in Fig ||, ||, and |J. 
Scenario A: The Cocoon at the Cluster Center. The 
pre-merger cluster is assumed to have an electron density 
of n c = 0.3 • 10~ 3 cm -3 and a temperature of kT = 3 keV 
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Fig. 2. Radio spectrum of the radio cocoon in scenario A 
at the end of phases 0-4. 

at the location of the radio galaxy. Due to the high en- 
vironmental pressure, we assume the internal pressure of 
the cocoon to be only twice the external pressure, and use 
Eq. [24| to obtain the initial volume Vq. Due to the strong 
synchrotron energy losses in this scenario, phase 2 can 
not last much longer than At2 = 0.1 Gyr, otherwise the 
revived fossil radio cocoon would not emit within the ob- 
servable radio frequency range. We assume that the shock 
wave of a cluster merger event increases the internal pres- 
sure by Pi, I Pi = 12 during the phase 3. This corresponds 
to a moderate shock with shock compression factor of 2.8, 
whereas strong non-relativistic shocks can have a com- 
pression factor of 4. A moderate shock is expected, since 
both merging clusters are expected to have temperatures 
of several keV and therefore sound velocities comparable 
to the merger velocity. 

As can be seen in Fig. || the compression caused by the 
merger shock wave gives rise to a burst of low frequency 
emission, but practically no high frequency emission. This 
is due to the rapid decay of the upper end of the elec- 
tron spectrum during phase 3, which essentially wipes out 
the adiabatic energy gains of these electrons. The source 
decays on a time-scale of a few tens of Myr, mostly due 
to the heavy synchrotron losses. If the radio cocoon is lo- 
cated in a more peripheral region of the cluster, where 
the density, the pressure and therefore the magnetic field 
strength inside the cocoon is much lower, these losses are 
also much milder. This lengthens the time scale over which 
the radiatively cooling synchrotron plasma can still be re- 
vived by the next passing shock, and thus rendered radio 
detectable. We, therefore, expect the radio relic phenom- 
ena to be found preferentially at larger cluster radii, and 



Fig. 3. Radio spectrum of the radio cocoon in scenario 
B at the end of phases 0-4. The luminosity at the end of 
phase 2 is too small in order to be displayed in this figure. 



less often near the cluster center (although projection can 
help some relics to appear near the cluster core) . The best 
environment to find cluster radio relics is, therefore, near 
the edges of the clusters. 

Scenario B: The Cocoon at the Cluster Boundary. 

The radio cocoon is assumed here to be born outside the 
cluster, in an environment of a dense galaxy filament, or a 
group of galaxies, say, with n c — 0.3 • 10 -5 cm -3 and kT = 
0.3 keV. The freshly injected radio plasma might be over- 
pressured by a factor of 100, leading to a short expansion 
phase. After this, the electrons within the expanded Mpc 
sized cocoon suffer mostly the IC-losses, allowing revival 
of the radio plasma even At 2 = 1 Gyr later. This can 
happen when the cocoon along with the ambient medium 
is crossed by the accretion shock of a cluster of galaxies, 
which might entail a pressure jump as large as P3/P2 = 
100, in order to heat the infalling cool gas to the cluster 
virial temperature of up to lOkeV. 

Scenario B can explain the steep and bent radio spec- 
trum of the cluster radio relic 0038-096 in Abell 85. 
An eye-fit to the radio spectrum (Fig. [3]) shows that 
the maximal electron momentum in this case is = 
10 4 (-B/VG) -1 / 2 . The magnetic field strength of the clus- 
ter relic was estimated from the minimum energy argu- 
ment to be B 



1 /xG (Feretti and Giovannini, 1996) 
and from the detection of excess X-ray emission at the 
location of the relic, which implies a field strength of 
B = 0.95 ± 0.10 ^G (|Bagchi et al., 1998|) if this emission 



refers to the IC scattered cosmic microwave background 
photons, otherwise a higher field strength. Using B — 1 /iG 
and p» = 10 4 and assuming a uniform environment with- 
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Fig. 4. Radio spectrum of the radio cocoon in scenario C 
at the end of phases 0-4. 



out expansion and compression, an age of 0.2 Gyr would 
result (Komissarov and Gubanov, 1994, and see Eq. 



0.001 




1000 



v/MHz 



But scenario B demonstrates that the radio plasma can be 
as old as 2 Gyr. This resolves the problem of the appar- 
ent cooling time of the electrons being too short for any 
nearby galaxy to have ejected the plasma and then moved 
to its present location with a typical velocity of a cluster 
member. For the long duration of phase 2 the resulting 
spectrum is fairly steep in the observable radio range. But 
this need not to be the case for a scenario with a shorter 
fossil phase. 

Srpnarin C: Thp Smoking Clin. Tn order to snhstan- 



tiate the last statement, we choose a. set of parameters 



Fig. 5. Radio spectrum of the radio relic in A85. The 
data w as compiled from the literature by Bagchi et al. 
( 1998 ). The solid line is a radio spectrum resulting from 
the canonical form (Eq. |l8| ) with a c = 3 and = 
11000 (B/^G)- 1 / 2 , and the dotted line one with a c = 2.65 
and V * = 10000 (S/^G)- 1 / 2 . 



pression at the cluster accretion shock wave might have 
taken a few hundreds of Myr. We assume a pressure jump 
of only P3/P2 = 40 at the shock wave, not higher, in or- 
der to allow the temperature of the post-shock gas t o stay 
below the a verage cluster temperature of 8.2 keV (Brie] 
ct al., 1992 ) . As can be seen in Fig. f|, the radio spectrum 



for scenario C which produces a cluster radio relic with 
relative flat, nearly unbent radio spectrum, like the relic 
1253+275 at the boundary of the Coma cluster of galaxies 
(Giovannini ct al., 1991, and see Fig. [l]). There, the prob- 
able source of the radio plasma is visible: it is the galaxy 
NGC 4789, which is located upstream of the relic, and 
which seems to be on an ascending orbit after a cluster 
core passage ( Enfilin ct al., 1998| ). Its narrow angle tailed 
radio outflow - the smoke of the gun NGC 4 789 - seems to 
be bent and dragged by the infalling matter (falling into 
Coma's gravitational potential) to the location of the relic. 
There, it brightens up and exhibits a steep but straight 
spectrum with a slope of 1.18. For details of the spectrum 



and th e geo metry, see Giovannini et al. (1991) and Enfilin 
et al. ( |1998|) . 

We follow a blob of the radio plasma, which might have 
been injected with an overpressure of a factor 10 into the 
infalling gas stream with possibly n c = 0.3 • 10 -5 cm~ 3 
and kT = 0.6 keV. The inflow of the plasma and the corn- 



below 1 GHz stays practically unbent for a couple of tens 
of Myr after the shock passage. 



7. Discussion 

We have argued here that adiabatic compression in 
cluster shocks can revive fossil radio plasma to radio 
detection, even up to 2 Gyr after its release from the 
parent radio galaxy. The computed radio spectra provide 
a good match to the observed spectra of cluster radio 
relics. Hence, this is a promising model for the regions of 
diffuse radio emission found in clusters of galaxies. Below, 
we summarize some merits and potential shortcoming of 
this model: 

Pros: 

— The observed connection of cluster radio relics to shock 
waves arises naturally in this model. 
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— The presence of the radio galaxy NGC 4789 and the 
morphological connection of its radio tails to the relic 
1253+275 seems to be the smoking gun of cluster radio 
relic formation by compression of fossil radio plasma. 

— The expected very high sound speed within radio 
plasma should virtually forbid shock waves of the am- 
bient medium to penetrate into radio cocoons. This 
renders adiabatic compression as a more plausible 
means of reviving the fossil plasma. 

— Cluster radio relics are rare. Not only a shock wave and 
fossil radio plasma have to be present in this model, in 
order to produce a cluster radio relic, but also the fossil 
radio plasma can not grow too old (of the order of 0.1 
Gyr in the center of clusters and 1 Gyr at their bound- 
aries). Otherwise the high energy end of the electron 
population would suffer too severe a depletion for the 
adiabatic gains during the compression to be able to 
shift them to radio emitting energi es. Since fossil radio 
plasma should be very common (Enfilin and Kaiser 
200^), this helps to explain the observed rarity of clus- 
ter radio relics. 

— The observed tendency of relics to appear at periph- 
eral positions in clusters follows in this model from the 
much shorter lifetimes of the radio emission and, conse- 
quently, much shorter ages which the revive-able fossil 
radio cocoons can have in the center of clusters, owing 
to the stronger synchrotron losses occurring there. 

Contras: 

— If the radio plasma has a high mass load, due to un- 
detectable cool gas, it should get shocked, ensuing 
Fermi-I shock acceleration to produce the observed ra- 
dio emitting electron population ( Enfilin et al., 1998| ). 



Still, the adiabatic gains of the relativistic particles are 
inescapable, although restricted in this case due to the 
limited compression factor of a shock wave. 
In the case when the fossil radio plasma is supported by 
a very hot, non-relativistic gas component, both adia- 
batic and shock acceleration are ineffective. The first 
due to the lower compressibility of a non-relativistic 
equation of state, while the second due to the inabil- 
ity of environmental shock waves to penetrate into the 
plasma. The fact that we do observe cluster radio relics 
tells us that this situation does not necessarily arise 
(or, alternatively, cluster radio relics are not related to 
radio ghosts). 

The giant cluster radio relic 2006-56 in Abell 3667 
shows a rim of flat radio emission (spectral index 
around 0.5) at its sharp outer edge (Rottgering et al. 
1990). This could be understood in a Fermi-I shock ac- 
celeration model due to the different ages of the elec- 
trons found at different distances from the energizing 



shock wave (Enfilin et al., 1998). There are two pos- 



sibilities to explain the observation of a flat spectrum 
rim within the adiabatic compression model: (a) there 
is a high mass load in the radio plasma, allowing the 



shock wave to successively compress the plasma, which 
should also produce a spatially different electron pop- 
ulation, and (b) the apparent rim of a flat spectral 
index is an artifact arising from an unmatched cover- 
age of the Fourier space in the radio interferometric 
observations at the two frequencies used. 

We conclude that the model presented here provides a 
fairly natural and promising explanation for the phe- 
nomenon of cluster radio relics. This model predicts the 
existence of a population of diffuse, ultra-steep, very low 
frequency radio sources inside and possibly also outside 
of clusters of galaxies, due to the age dependence of the 
upper frequency cutoff of radio emission arising from the 
revivable fossil radio plasma. The ongoing rapid improve- 
ments in the sensitivities at low radio frequencies (Giant 
Meterwave Radio Telescope: (GMRT), Low Frequency Ar- 
ray (LOFAR)) may therefore open a new window on cos- 
mological structure formation by detecting shock waves 
marked by such relic radio sources. 
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